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1 Introduction 

The Backhand transformation equations play an important role in the study of integrable equations JjJ . 
They generate the Backlund transformation which allows one to obtain the exact N soliton solution of the 
nonlinear equation, and also provide the Lax pair. The integrable equation itself is given as a consistency 
condition of the Backlund transformation equations in these schemes. 

On the one hand, the cellular automaton models which have exact N soliton solution were found in 
the last decade JsJ . In many cases, these soliton cellular automaton models are obtained as a special limit, 
known as the ultradiscrete limit, of the discrete integrable equations. Typically, the limiting procedure 
simply replace the summation by the Max operator and the product by the summation. Despite of this 
simple correspondence, this procedure causes some ambiguities concerning to the Max operator. For 
example, we can remove the variable x from the equation a + x = b + x, however, we can not remove the 
variable x from 

Max(a,x) = Max(6, x). 

These type of difficulties prevent us to consider the consistency condition of the Backlund transformation 
equations for the soliton cellular automaton model. 

In this paper, we consider the ultradiscrete KP equation and its Backlund transformation equations, 
since it reduces to many soliton cellular automaton models. The ultradiscrete KP equation was first 
proposed in the paper |^|, as the unified equation of a class of generalized Box and Ball systems. It also 
contains the ultradiscrete Toda equation etc. 

We show the consistency condition of the Backlund transformation equations of the ultradiscrete 
KP equation as follows. First, we propose an extension of the Backlund transformation equations of 
the discrete KP equation An important aspect of the extended Backlund transformation equations 
is that, it can be expressed as the vanishing of the product of a matrix and a vector. The consistency 
condition becomes the vanishing of the determinant of the matrix. This will be shown in the next section. 
Second, in section 3, we show the algorithm to eliminate N variables from N ultradiscrete linear equations, 
before considering the ultradiscretization of the discrete KP equation. Finally, in section 4, we derive 
the Backlund transformation equations of the ultradiscrete KP equation, starting from the extended 
Backlund transformation equations of section 2. By applying the algorithm of section 3, we derive the 
consistency condition of the Backlund transformation equations. Some examples of soliton solutions are 
also shown. 
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2 Backlund transformation equation for the discrete KP equa- 
tion 



Let us begin with writing down the discrete KP equation for a function / of three discrete variables p, q, r 
@ 0, 

ZpZqr fpfqr ZqZrpfqfpr ~t~ %r Zpqfr fpq — 0- (1) 

Here, ZiS are arbitrary complex constants and Zij = Zi — Zj. The lower subscripts of / indicate to increase 
the corresponding variables by one, 

f P = f{p+i,q, r), 

f pq = f(p+l,q+l,r). (2) 

Despite of the simplicity this equation plays an important role in the study of the integrable equations. 
It reduces to many integrable equations known in the physics and mathematics and is satisfied by 
the transfer matrices of a class of solvable lattice models || , || , @ , and is also satisfied by the amplitude 
of a string model 0] . 

The Backlund transformation equations for the discrete KP equation are a system of linear equations 

Zrfr9q z qfq9r "i" ZqrfqrfJ 0> 

Zrfrdp Zpfpfjr ~t~ Zp r fp r g — 0. (3) 

If we substitute a solution of the discrete KP equation into the function / in this equation, then we can 
show that there exist nontrivial solutions g, which satisfy the discrete KP equation. The compatibility 
condition of these two equations is in the central of this scheme. It guarantees the existence of the solution 
for these equations, and it also guarantees that the solution, g also satisfies the discrete KP equation. 
See jl0| and || for example. However, the consistency condition does not work for the ultradiscretized 
system, because of the ambiguity of the Max operator stated above. 

In this paper, we extend the Backlund transformation equations (Q), by considering the additional 
bilinear equations deduced from the discrete KP equation and its Backlund transformation equations (|J) . 
The additional bilinear equations provide another expression of the Backlund transformation equations. 

Let us rewrite the Backlund transformation equations, by dividing them by g and move the terms 
other than f qr and f pr to the right hand side 

~Z r f r gq + Zqfgg r 



fqr 



Zq r g 



r z rfrgp + z pfp9r , .\ 

Jpr — ■ 
^p r y 
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Substituting these expressions of f qr and f pr into the discrete KP equation, we can obtain an additional 
bilinear equation 

z pfp9q + z qfq9p ~ z pqfpq9 = 0- (5) 

Similarly, substituting f p and f q appearing in the equations ([|) into the discrete KP equation, we obtain 
another bilinear equation 

z pqfpq9r ~t~ z qr fqr9p z pr fpr9q — 0. (6) 

Notice that the equations (^), (0) and (^) are linearly independent, unless the function / satisfies the 
discrete KP equation. These four equations are the extension of the Backlund transformation equations, 
which we employ in this paper. 

An important aspect of the extended Backlund transformation equations is that they can be summa- 
rized into the product of a matrix and a vector. Indeed, the extended Backlund transformation equations 
consist of four equations and the number of g appearing in the same equations (i.e g, g p , g q , g r ) is four, 
thus, we can write down the equations by using a 4 x 4 matrix and a 4 component vector 



z r fr z qfq z rqfrq 

Z r f r z pfp z pr fpr 

Z qfq Z pfp z qpfpq 



Zqrfqr z rpfpr z pqfp 









9 P 




9q 




9r 
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0. 



(7) 



The consistency condition for these equations is the existence of the nontrivial solution, g and it is nothing 
but the vanishing of the determinant of this matrix. In this case, the matrix is antisymmetric matrix, 
thus, we can use the fact that the determinant of the antisymmetric matrix becomes the square of a 
pfaffian. The vanishing of the pfaffian becomes just the discrete KP equation. 

Since the consistency condition is very different from the consistency condition for the original 
Backlund transformation equations, there is a chance to avoid the difficulty of the Max operator. 

3 Ultradiscrete linear equations 

Of course, we can remove the Max operator from the ultradiscrete Backlund transformation equations, 
by dividing it into cases. However, it is complicated and requires the individual consideration. Instead 
to do so, we study the generic consistency condition for the ultradiscrete linear equations. We show an 
algorithm to eliminate variables, starting from the consideration of two linear equations. 

Before proceeding to the generic cases, let us consider the following two examples of two ultradiscrete 
linear equations. Our first example is the following couple of equations 



a + x = Max(a + x, b + y), 
c + x = Max(c + x, d + y). 



(8) 



Here, a, 6, c, d, a, c are the constant coefficients, and x, ?/ are the variables of these two equations. We can 
eliminate these variables, by considering the trivial identity 

Max(Max(a + x,b + y) + d, c + b + x) = Max(a + d + x, Max(c + x,d + y) + b). (9) 

Substituting the first and second equations of (||) into the left hand side and the right hand side of this 
equation respectively, we obtain the relation 

M&x(d + d, c + b) = Max(a + d, c + b). (10) 

This relation should be satisfied, if there are finite solutions, x and y. Thus, we can eliminate all variables, 
without dividing Max operators into the cases. 

Our second example is the following couple of linear equations 

a + x — Max(a + x, 6 + y), 

d + y = Max(c + x,d + y). (11) 

We can also eliminate all variables from these equations without dividing them into cases. Indeed, adding 
these two equations and applying these two equations to the left hand side once again, we obtain the 
relations 

a + d + x + y 

= Max(a + c + 2x, a + d + x + y, b + c + x + y, b + d + 2y) 

= Max(a + d + x + y,a + d+x + y,b + c + x + y). (12) 

We can eliminate all variables, by subtracting the x + y from the first part and the third part of these 
equations, 

Max(a + d, a + d, b + c) = a + d. (13) 

In the general two linear equations cases, we can not eliminate all variables without dividing them into 
cases. However, previous two examples make the consideration simple. The most general ultradiscrete 
two linear equations for two variables x and y are written as 

Max(a + x, b + y) — Max(a + x, b + y), 

Max(c + x,d + y) = Max(c + x.d + y). (14) 

For example, if a + x in the left hand side of the first equation is larger than b + y and c + x in second 
equation is larger than d + y, these equations become as 

a + x = Max(a + x, b + y), 

c + x = Max(c + x,d + y). (15) 
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This is the case of our first example. By applying the result of the first example, we can eliminate all 
variables from these equations. 

Max(a + d, c + b) = Max(a + d, c + b) (16) 

Applying the similar argument to the rest of the three cases 

a + x > b + y and c + x > d + y =>• Max(a + d, c + b) — Max(a + d, c + b) , 

a + x < b + y and c + x < d + y =>• Max(6 + c,d + a) = Max(6 + c, d + a) , 

a + x > b + y and c + x < d + y =>■ d + d = Max(a + d, a + d, b + c), 

a + x < b + y and c + x > d + y =3- b + c = Max(6+ c, b + c, a + d) , (17) 

and summarizing these results in one equation, we can obtain the relation 

Max(a + d, a + d, b + c, b + c) = Max(a + d, a + d, b + c, b + c). (18) 

One can easily confirm that this relation is satisfied in all four cases. This relation should be satisfied, if 
there are finite solution, x and y for the generic two linear equations (|l4|). 
On the one hand, we can rewrite the equations ( |l4| ) as 

Max(a + z,b) = Max(a + z,b), 

Max(c + z,d) = Max(c + z,d). (19) 



Here, we subtract y from equations (14), and replace x — y by z. Apparently, if there is finite z which 
satisfies these equations, the relation ( |l8| ) should be satisfied. Since, a, 6, c, d, a, 6, c, d are arbitrary con- 
stants, we can eliminate one variable from two linear equations of any number of variables. Applying 
this reduction recursively, we can reduce N linear equations for N variables to 2 linear equations for 
2 variables, and the rest of 2 variables can be eliminated from two linear equations. So that, we can 
eliminate N variables from N linear equations. 

In the next section, we consider the consistency condition of the Backlund transformation equations 
for the ultradiscrete KP equation by using this algorithm. 

4 The Backlund transformation equations for the ultradiscrete 
KP equation 

Let us write the discrete KP equation, once again. It was the following equation 

ZpZq r fpfq r 4~ Z r Zpqf r fpq — ZqZp r fqfpr ■ V^^O 



Here, / is a complex valued function in the original discrete KP equation, but we restrict the value to 
positive real numbers, to take the ultradiscrete limit. We also restrict the values of z^s to real numbers, 
and assume z p > z q > z r since it does not lose the generality. Notice that in this choice of coefficients, 
both of the left hand side and the right hand side of the equation ( po|) become positive. 

The ultradiscrete limit of the discrete KP equation is taken as follows. First, we transform the 
dependent variable / and the coefficients Zi, into F and Zi, Z^j, as 

/ = exp(— ) , z l = exp(^-) , Zij=exp(— ) . (21) 



Here, e is a positive real parameter. Second, we apply the log function to both side of the equation (|20|) 
after multiplication of e, and take the e — > +0 limit. By using the following two identities, 

A B 
lim loge(exp( — ) + exp( — )) = Max(A B), 
e-»+o e e 

A B 

lim loge(exp(— ) x cxp( — )) = A + B, (22) 
e^+o e e 

we can obtain the ultradiscrete KP equation from equation J20|). 

Z q + Z pr + F q + F pr = Max(.Zp + Z qr + F p + F qr , Z r + Z pq + F r + F pq ) (23) 

As in the case of the discrete KP equation, this equation reduces to many soliton cellular automaton 
models. See ||, for example. 

The Backlund transformation equations for the ultradiscrete KP equation are obtained by applying 
the similar procedure to the extended Backlund transformation equations. Applying the transformation 
@ and 

<7 = exp(-), (24) 

and taking the e — > +0 limit, we obtain the equations 

Z q + F q + G r = Max[Z r +F r + G q , Z qr + F qr + G] , 
Z p + F p + G r = Max[Z r + F r + G p , Z pr + F pr + G], 
Zp + F p + Gq = Max[Z 9 + F q + G p , Z pq + F pq + G], 
Z pr + Fp 7 . + G q = Max[Z qr + Fqr + G p , Z pq + F pq + G r ], (25) 

The consistency condition for these equations is obtained by eliminating the all variables from these 
equations. Let us notice that we can eliminate G r and G g , by substituting second equation and third 
equation to first equation and fourth equation, respectively. Then, we obtain the following two linear 
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equations for two variables G p and G 

Max(Z r + Z q + F r + F q + G p ,Z q + Z pr + F q + F pr + G) 

= Max(Z r + Z q + F r + F q + G p , Max(Z r + Z pq + F r + F pq , Z qr + Z p + F qr + F p ) + G)), 

Max(Z pr + Z q + F pr + F q + G p , Zpr + Z pq + Fp 

r H~ ^pq 

+ G) 

= Max(Max(Z 9r + Z p + F qr + F p , Z pq + Z r + F pq + F r ) + G p , Z pq + Z pr + F pq + F pr + G). 

(26) 

By applying the relation (^), we can eliminate the rest of two variables. After some calculation, we can 
obtain the equation for F 

Max(Z q + Z r + Z pq + Z pr + F q + F r + F pq + Fpr, 

Zp r + Z q + Fp r + F q + Max(Z r + Zp q + F r + Fpq, Z qr + Z p + Fqr + F p ), 

Z q + Zpr + F q + F pr + Max(Z gr + Zp + Fqr + F p , Z pq + Z r + F pq + F r )) 

= Mdx(Z q + Z r + Zpq + Z pr + F q + F r + F pq + F pr , 

2{Z pr + Z q + Fpr + Fq), 2Max(Z r + Z pq + F r + F pq , Z qr + Z p + F qr + F p )). (27) 

Notice that we can remove the first terms of the both hand side, since they are smaller than the other 
terms. Moving the rest of the right hand side to the left hand side, we obtain the equation 

\\Z q + Zpr +F q + F pr ~ Max(Z p + Zqr + F p + F qr , Z r + Z pq +F r + Fpq\\ = 0. (28) 

This is nothing but the absolute value of the ultradiscrete KP equation. 

Since we can change the roles of F and G, by shifting each of the equations of ( p5|) appropriately, G 
should also satisfy the ultradiscrete KP equation. Thus, we can show that equations ( p5|) really generate 
the Backlund transformation of the ultradiscrete KP equation. 

In the rest of this section, we derive the two soliton solution and conjecture the form of the N soliton 
solution. Let us assume the existence of the trivial solution, F = G = 0, then the conditions for the 
coefficients of equations (|25|), 

Zq = Max[Z r , Zqr], 
Zp = Max[iv r , Zqr, Zpq], 

Zpr=M&x[Zqr,Zpq], (29) 

should be satisfied. Notice that Z r , Z qr and Z pq determine the rest of the coefficients. We rewrite these 
three coefficients as 

Zpq = Z\ , Zq r — Z2 , Z r = Z3 , (30) 
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for simplicity. Furthermore we assume the relations Z\ < Z 2 < Z%, then the rest of the coefficients 
becomes as 

Z q — Z 2 , Z p = Z\ , Z pr — Z\. (31) 

Of course we can consider the other relations for Z\, Z 2 and Z3, however, it does not change the following 
argument. 

The one soliton solution is derived as follows. First, we substitute the trivial solution F = into the 
Backlund transformation equations (p5[) , and assume the form of G as 



G = a ■ p. (32) 
Here, a is the three dimensional vector of constant coefficients, and p is the vector of coordinates p, q, r 

a= (a p ,a g ,a r ) , p=(p,q,r). (33) 



The dot between them represents the vector product. By substituting the G of equation (|32|) into the 
Backlund transformation equations (^5|), we obtain the following dispersion relations 

Z 2 + a r = Max[Z 3 + a q , Z%] , 
Z\ + a r = Max[Z 3 + a p , Z\], 
Zi+ a q = Max\Z 2 + a p , Z\\, 
Zi + a q = Max[Z 2 + a pi Zi + a r }. (34) 

We can solve these dispersion relations as 

a p = a + Z\ , oi q — Max(a + Z 2 , 0) , a r — Max(a + Z3, 0). (35) 

Here, a is an arbitrary real parameter. The one soliton solution of the ultradiscrete KP equation becomes 
the superposition of these representations of G which is written in "Max" operator, 

G = Max(a -p , (3-p). (36) 

where (3 is a three dimensional vector ((3 p , (3 qi (3 r ) which satisfies the dispersion relations. 

Figure 1 shows the time evolution of three soliton solution. We transform the dependent variable / 
and coordinates p, q, r into U and l,m,n as 

U (p, q, r) = f(p, q,r + l)+f(p + l,q+ 1, r) -f(p+l,q,r)- f(p, q + 1, r + 1), 

I = p — r , m = q , 11 = p + r. (37) 



The / axis is parallel to the horizontal line and m axis is parallel to the vertical line. 

We see the propagation of one soliton solution in the large I and large m region, in Fig. 1. 

Similarly, we can obtain the two soliton solution by substituting the one soliton solution ( |36| ) to the 
function F in the Backlund transformation equations, and solve the equations for G. Let us assume the 
form of G as 

G = Max(c Q + a ■ p + 7 • p , cp + (3 ■ p + 7 • p), (38) 

where c a and cp are the constants which should be determined later, and 7 is also a three dimensional vec- 
tor which satisfy the dispersion relations. Substituting this expression into the Backlund transformation 
equations, we obtain the following relations of c a and cp 

c a = cp, if 7<«. P, 

c a - a = cp - (3, if a , (3 < 7, (39) 

There is no solution of c a and cp if a < 7 < (3 or (3 < 7 < a. 

These results can be summarized into the following forms, without the loss of the generality 

c a = Max(a, 7), 

cp = Max(a, 7) , (a , (3 < 7 or 7 < a , (3). (40) 

The two soliton solution of the ultradiscrete KP equation becomes the superposition of these solutions, 

G = Max(Max(a, 7) + a ■ p + 7 • p , Max(a, S) + a ■ p + 5 ■ p , 

Max(/3,7) +0-p + j-p , Max(/3,7) +0-p + 5-p) , (41) 

where 7 and S should satisfy one of the following conditions 

7 , <5 < a , (3 
7 <a , 13 < 8 

a, 13 < l,S. (42) 

These arguments allow one to suppose the form of the N soliton solution. The N soliton solution may 
have the following form 

/ = Max w (J2 Max« j , a? ) + £ dT • £)■ (43) 

i<j i 

Here, {a} = (a\, 02, un) and oVs take value of 1 or 0. Max^} means to take the maximal value of all 
the combination of the cr's. The vectors d*j° and cti are constant coefficients which satisfy the dispersion 
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relations, and we also assume that one of the following conditions is satisfied for all combinations of i < j 

a° , a- < a° , a], 
a- < oP 3 , a] < a], 
a° , a] <a° , a*. (44) 

We do not have the proof of the expression (|43|), but we confirm that it really satisfies the Backlund 
transformation equations till N = 4. The complete proof will be given in the other paper. 

5 Conclusion 

We consider the consistency condition for the Backlund transformation equations of the ultradiscrete 
KP equation. We extend the Backlund transformation equations for the discrete KP equation into the 
form of the product of a matrix and a vector. The consistency condition becomes the vanishing of the 
determinant of the matrix. In this case, the matrix becomes an antisymmetric matrix, and we can obtain 
the discrete KP equation by using the relation between the determinant of an anti symmetric matrix and 
a pfaffian. We also show the algorithm to eliminate the variables from the ultradiscrete linear equations. 
By using this algorithm, we obtain the consistency condition of the Backlund transformation equations of 
the ultradiscrete KP equation. The consistency condition becomes the absolute value of the ultradiscrete 
KP equation, instead of the square of the pfaffian in the discrete KP equation case. 
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Figure 1: Time evolution of 3 soliton solution 
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